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Abstract: In sharp contrast to the corresponding clas- 
sical systems cases it is not yet understood how to define 
a mechanical quantity with the interpretation of entropy 
creation rate for nonequilibrum stationary states of finite 
quantum systems with finite thermostats. Some aspects 
of this problem are discussed here in cases in which iden- 
tifying entropy creation rate as a mechanical observable 
might be possible. 
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1. Classical systems and thermostats 

The aim of this paper is to propose a notion of ther- 
mostat acting reversibly on a quantum mechanical sys- 
tem which is analogous to the Gaussian or Nose-Hoover 
thermostats that have been so important in recent times 
for the development of research non nonequilibrium. At 
the same time we remark that "entropy creation" can be 
associated with the thermostats discussed here and can 
be given a meaning that makes it susceptible of experi- 
mental measurements. This is important also because it 
allow us to extend the "chaotic hypothesis" to such sys- 
tems and to test the "fluctuation relation" that holds, at 
least as a formal consequence. I begin by recalling the no- 
tion of classical thermostat that will be genreralized here 
and some results concerning classical nonequilibrium sta- 
tistical mechanics. 

A thermostat is, physically, a device that extracts heat 
created inside a mechanical system subject to non conser- 
vative forces, thus allowing control of the energy build- 
up: in this way a forced system can reach a stationary 
state. Such a state is however deeply different from an 
equilibrium state. Heat, matter or electric currents may 
be present and the dissipation associated with the ther- 
mostats implies that the probability distribution \x that 
describes the statistical properties of the system is not 
even close to the familiar Maxwell-Boltzmann distribu- 
tions for equilibrium. This feature makes the problem 
of the statistics of stationary noncquilibria particularly 
interesting and challenging. 

In the classical case a rather general model for a system 
in contact with thermostats is represented in Fig. 1 , [1] . 
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Fig.lReservoirs occupy finite regions outside Co, e.g. sectors 
Cj C R , j = 1,2.... Their particles (with mass 1) are con- 
strained to have a total kinetic energy Kj constant, by suitable 
forces, so that the reservoirs "temperatures" Ti are well defined. 
Fixed j the label i is in 1, . . . ,Ns. 

The equations of motion for the particles positions 
X ,Xj, . . . , X n are written in Fig.l in terms of 

(1) The potential energies for the Nj particles inside each 
region Cj, Uj (X ■ ) , which will be assumed bounded, for 
simplicity; all masses are m = 1, also for simplicity, 

(2) The potential energies between particles in Cq and in 
Cf Uj(Xn. Xj). also assumed bounded, 

(3) The external, non conservative, positional forces F 
acting on the particles in Co, 

(4) The thermostat forces o.jX i which are so defined that 

the total kinetic energy Kj = £Xj in each thermostat is 
strictly constant. Such forces will be imagined realized 
by imposing constancy of Kj via Gauss' principle, [2, 
Sec.9]. This gives, [1], 
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where Wj is the work done per unit time by the system 
particles on the j— th thermostat particles and Uj is the 
internal potential energy of the j-th reservoir: Wj = 
—Xj ■ dx Uj(X q, X j). Thus Wj will be identified with 
the heat Qj ceded per unit time by the system to the j-th 
thermostat. 

(5) The value of the constant Kj will be written Kj — 
ikNjksTj {ks = Boltzmann's constant) and will define 
the temperature of the j-th thermostat, [1]. 

A brief computation yields the divergence of the equa- 
tions of motion, i. e. of the phase space contraction rate, 
for the velocity-position coordinates X_, X_ 
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where, remarkably, e(X, X] can be interpreted as the 
entropy creation rate 
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Eq.(1.3) are correct up to 0(N : ) if N = miniVj as the 
addends should contain also a factor (1 — ^-) to be exact: 
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for simplicity 0(l/N) corrections will be ignored (their 
inclusion would imply trivial changes without affecting 
the physical interpretation), [1]. 

An additive total derivative, R(2Q in this case, of a 
bounded quantity does not affect the long time fluctua- 
tions. Therefore the average phase space contraction and 
the average entropy creation rate have the same average 
a + = e + and, assuming e + ^ 0, the same large devi- 
ations rate function ((p) for p = -^-p J Q T cr(S t (X_, X_))dt 

and for ^ ^ e{S t {X,X))dt. 

Furthermore the equations of motion are reversible so 
that, under the Chaotic Hypothesis, the Fluctuation The- 
orem yields, [f, 3-6], the symmetry property 



C(-P) = C(p) ~ P£+, \p\< P* 
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for the large deviations rate ((p) of p. The identi- 
fication between phase space contraction and entropy 
creation rate is thus motivated, ft should be noted 
that Eq.(1.3) shows that even in experiments, a case in 
which one hardly knows the equations of motion and the 
phase space divergence, the divergence can be measured 
through heat flow measurements, at least as far as its fluc- 
tuations over large time are concerned, and give nontriv- 
ial consequences like the fluctuation relation, Eq.(1.4). 

2. Quantum systems and thermostats 

Is it possible to formulate a dissipation theory analo- 
gous to the one developed for classical systems when the 
quantum nature of the system in Co cannot be neglected? 

At first it might seem almost impossible: in quantum 
systems average kinetic energy is not identified with tem- 
perature; and all motions are quasi periodic if the sys- 
tem is of finite size (as in our examples in Sec.l), so that 
strictly speaking no chaos is possible. 

A way out, explored in the literature, would be to 
imagine the thermostats as infinite systems whose state 
far from Co is a Gibbs state at a well defined temperature, 
[7-10]. This is a point of view that could also be taken in 
the classical case: however the recent progress in classical 
statistical mechanics was sparked by the introduction of 
finite size thermostat models and this is the path that 
will be attempted here. 

Thermostats have, usually, a macroscopic phenomeno- 
logical nature: in a way they should be regarded as clas- 
sical macroscopic objects. Therefore it seems natural to 
model them as such: thus their temperature can be de- 
fined as the average kinetic energy and the question of 
how to define it does not arise. 

Consider the system in Fig.l when the quantum na- 
ture of the particles in Co cannot be neglected. Suppose 
first that the nonconservative force FQCq) acting on the 
system vanishes, i.e. consider the problem of heat flow 
through Co- Let H be the operator on L 2 (Cg Ar °), space 
of symmetric or antisymmetric wave functions \E', 
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j>0 

where Ax is the Laplacian, and note that its spectrum 
consists of eigenvalues E n — E n ({X_j}j >a ), for 2Lj fixed. 

A system-reservoirs model can be the dynamical sys- 
tem on the space of the variables ({2Lj}, {2Lj})j>o) 
defined by the equations (where (-)^ is the expectation 
in the state 'J') 



-ift*C£o) = (#*)QEo). and for 3 > 
Zj = -(djUiQCj) + (djUj^JQ))^ aj Xj (2.2) 
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here the first equation is Schrodingcr's equation, the sec- 
ond is an equation of motion for the thermostats parti- 
cles, [1], similar to the one in Fig.l, whose notation for 
the particles labels is adopted here too. Evolution main- 

tains the thermostats kinetic energies Kj = ^X_^ exactly 
constant so that they can be used to define the ther- 
mostats temperatures Tj via Kj — ^ksTjNj, as in the 
classical case. 

Let jio({d^}) be the formal measure on L>2{C iNo ) 
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with <l/ r , real and imaginary parts of ^. The formal 
phase space volume element /^({(M}) x v(dlfi) with 



v{dXdX) d = f J[ (dXj dXj (5(aJ - SNjkBTjj^J (2.4) 

j>0 

is conserved, by the unitary property of the wave func- 
tions evolution, just as in the classical case, up to the 
volume contraction in the thermostats, [1]. 

def 

If Qj = (Wj)^ and R is as in Eq.(1.3) the contrac- 
tion rate a of the volume element in Eq.(2.4) is given by 
Eq.(1.2) with e, that will be called entropy creation rate, 
defined by Eq.(1.3). 

In general solutions of Eq.(2.2) will not be quasi pe- 
riodic and the Chaotic Hypothesis, [2, 11], can be as- 
sumed: if so the dynamics should select an invariant 
distribution fi. The distribution will give the statis- 
tical properties of the stationary states reached starting 
the motion in a thermostat configuration (2Lj,2Lj)j>o 

randomly chosen with "uniform distribution" v on the 
• 2 

spheres X_^ = SNjIcbTj and in a random eigenstate of 
H. The distribution \i, if existing and unique, could be 
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named the SRB distribution corresponding to the chaotic 
motions of Eq.(2.2). 

In the case of a system interacting with a single ther- 
mostat the latter distribution should be equivalent to the 
canonical distribution. 

Hence an important consistency check for the model 
proposed in Eq.(2.2) is that there should exist at least 
one stationary distribution equivalent to the canonical 
distribution at the appropriate temperature T\ associ- 
ated with the (constant) kinetic energy of the thermo- 
stat: Ki = jk B TiNi. In the classical case this is an 
established result, [12], [1, 2]. 

Hence an important consistency check for the model 
proposed in Eq.(2.2) is that there should exist at least 
one stationary distribution \i equivalent to the canonical 
distribution at the appropriate temperature T\ associ- 
ated with the (constant) kinetic energy of the thermo- 
stat: K\ = ^ksTiNi. In the classical case this is an 
established result, [12], [1, 2]. 

The check should be performed also in the present case, 
thereby providing further motivation and support for the 
model in Eq.(2.2). A first candidate for /j, might be to 
attribute a probability proportional to dX_ x dX_ x times 

oo 

]T e"^»J(* - * n (2d) e*») d<p n 6(£ - 2Ki) (2.5) 

71=1 

where ^ is the wave functions for the system in C and 
X\ , X\ are positions and velocities of the thermostat par- 
ticles and ip n G [0, 2ir] is a phase, E n — E n (X_ 1 ) is 
the n-th level of H(X_ 1 ) with ^f n (K.i) the correspond- 
ing eigenfunction. The relation between the distribution 
in Eq.(2.5) and EIG(p) in [13] should be noted. However, 
as pointed out by a referee, Eq.(2.5) is not invariant un- 
der the evolution Eq.(2.2) and it seems difficult to exhibit 
explicitly an invariant distribution. 

3. Some consequences. Conclusion. 

The simplest case arises when V = 0: i.e. no noncon- 
servative force acts on C and Eq.(2.2) models heat flow, 
through Co, between the various reservoirs. 

Solutions of Eq.(2.2) are also reversible: time reversal 
being the change in sign of the velocities and the con- 
jugation of the wave function ^(Aq). Hence under the 
Chaotic Hypothesis the Fluctuation Theorem, [3], [2], see 
Eq.(1.4), would hold for the entropy creation rate fluctu- 
ations in the SRB distribution: however since the phase 



space is infinite dimensional corrections to Eq.(1.4) have 
to be expected for large p, as in [4], and should be dis- 
cussed on a case by case basis. Note that the fluctuation 
theorem extension is also here immediate once a model 
is properly formulated: as it was in the analogous cases 
of infinite thermostats, [8]. 

If a nonconservative force F(Xq) acts on the sys- 
tem and has a (multivalued) potential V(X ). so that 
E-iXjc) = — dyiXfi) is the force on the particles in Co, 
then in absence of thermostats generally the system will 
not reach a stationary state: this is true both in the clas- 
sical and in the quantum cases. In the latter case the 
Schrodinger equation, with H modified into H + V, will 
not have eigenvectors because of the multivaluedness of 
the potential V: actually it will not even be well de- 
fined. It should, however, be interpreted as an equation 
for the wave function defined on the "covering space" 
Q of C™ in which the potential V(X ) becomes sin- 
gle valued, and the thermostats should have the effect of 
allowing reaching a stationary state described by wave 
functions *„Q£o + £(£)) with G L 2 {Ct) and £(t) a 
suitable "flow". 

Whether this really happens is, however, an open prob- 
lem even in the classical case; there it has been called the 
problem of efficiency of the thermostats, [1, 14], and it 
has been studied only in a few numerical simulations in- 
volving long range or short range particles interactions. 
A fortiori it is an open problem in the quantum case. 

Finally it should be stressed that Eq.(2.2) provides a 
model of finite thermostat for a quantum system and 
therefore may be suited for simulations and tests. 

Identification of phase space contraction rate as the en- 
tropy creation rate (up to an additive total time deriva- 
tive) is an achievement of the recent research in classical 
nonequilibrium statistical mechanics, [12], which has led 
to the possibility of nonequilibrium simulations without 
the need of considering infinite thermostats and, subse- 
quently, to general results like the Fluctuation Theorems, 
[3], [15, 16], and the possibility to test them experimen- 
tally, [17], [1, 18], and even to make use of them in appli- 
cations. In this note the proposal that rather straightfor- 
ward extensions to quantum nonequilibrium are possible 
has been discussed. 

Acknowledgement: I thank F. Zamponi for pointing 
out reference [8]. 



[1] G. Gallavotti, Chaos 16, 043114 (+6) (2006). 

[2] G. Gallavotti, Statistical Mechanics. A short treatise 
(Springer Verlag, Berlin, 2000). 

[3] G. Gallavotti and E. G. D. Cohen, Physical Review Let- 
ters 74, 2694 (1995). 

[4] F. Bonetto, G. Gallavotti, A. Giuliani, and F. Zamponi, 



Journal of Statistical Physics 123, 39 (2006). 

[5] G. Gallavotti, Mathematical Physics Electronic Journal 
(MPEJ) 1, 1 (1995). 

[6] D. Ruelle, Journal of Statistical Physics 95, 393 (1999). 

[7] J. P. Eckmann, C. A. Pillet, and L. R. Bellet, Commu- 
nications in Mathematical Physics 201, 657 (1999). 



4 



[8] J. Kurchan, cond-mat/0007360 (2000). 
[9] D. Ruelle, Journal of Statistical Physics 98, 55 (2000). 
[10] D. Ruelle, Communications in Mathematical Physics 
224, 3 (2001). 

[11] G. Gallavotti and E. G. D. Cohen, Journal of Statistical 

Physics 80, 931 (1995). 
[12] D. J. Evans and G. P. Morriss, Statistical Mechanics of 

Noneqmhbrium Fluids (Academic Press, 1990). 
[13] S. Goldstein, J. Lebowitz, R. Tumulka, and N.Zanghi, 

Journal of Statistical Physics 125, 1193 (2006). 
[14] P. Garrido and G. Gallavotti, Journal of Statistical 

Physics cond-mat/0607231 and mp_arc 06-197, doi: 



10.1007/sl0955 (2007). 
[15] G. Gallavotti, Annales de 1' Institut H. Poincare 70, 429 

(1999 and chao-dyn/9703007). 
[16] E. G. D. Cohen and G. Gallavotti, Journal of Statistical 

Physics 96, 1343 (1999). 
[17] F. Bonetto, G. Gallavotti, A. Giuliani, and F. Zamponi, 

Journal of Statistical Mechanics p. P05009 (2006). 
[18] G. Gallavotti, Journal of Statistical Mechanics (JSTAT) 

pp. P10011 (+9) (2006). 

e-mail: giovanni . gallavotti@romal . inf n . it 

web: http : //ipparco . romal . inf n . it ReVtEx 



